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Abstract: The exact mass gap of the 0{N) Gross-Neveu model is known, for arbitrary 
N ^ from non-perturbative methods. However, a "naive" perturbative expansion of the pole 
mass exhibits an infinite set of infrared renormalons at order formally similar to the 
QCD heavy quark pole mass renormalons, potentially leading to large 0(A) perturbative 
ambiguities. We examine the precise vanishing mechanism of such infrared renormalons, 
which avoids this (only apparent) contradiction, and operates without need of (Borel) sum- 
mation contour prescription, usually preventing unambiguous separation of perturbative 
contributions. As a consequence we stress the direct Borel summability of the (l/N) per- 
turbative expansion of the mass gap. We briefly speculate on a possible similar behaviour 
of analogous non-perturbative QCD quantities. 
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1. Introduction 

The (1 + 1) dimensional 0{N) Gross-Neveu (GN) model[l, 2] often serves as a simpler toy 
model for more complicated theories like QCD, sharing with it the properties of asymptotic 
freedom and dynamical mass generation, while being an integrable model with many exact 
results available. The exact mass gap (associated with the breaking of the discrete chi- 
ral symmetry), for arbitrary N, has been evaluated from non-perturbative (NP) methods, 
more precisely from exact S matrix results [3] associated with Thermodynamic Bethe Ansatz 
(TEA) methods [4]. Yet from a more general viewpoint, it can be an interesting issue (since 
still not fully clarified, in our opinion) to study the precise matching between those two- 
dimensional NP exact results on one side and the standard perturbative behaviour on the 
other side. This may give some more insight on the short/long distance physics interplay 
for more involved theories like four-dimensional QCD. In particular one apparent puzzle 
arises, once realizing, as we examine here, that the (naive) perturbative expansion of the 
pole mass of the (massive) GN model suffers at next-to-leading 1/7V order from potential 
ambiguities, due to the presence of severe infrared renormalons, which are indeed formally 
completely similar to the quark pole mass renormalons [5]. While the TEA GN mass gap [4], 
and a fortiori its next-to-leading expansion[6], are unambiguously determined. Ac- 
tually, independently of TEA results, it is expected on general grounds that any truly 
NP calculation should be free of such ambiguities [7], i.e. that infrared renormalons are 
perturbative artifacts. Eut up to now only a few NP results have been explored from this 
perspective, even for integrable models, thus we find instructive to examine in some details 
how exactly the NP contributions, here fully controllable at least at order, organize 
themselves so that the necessary vanishing of such ambiguities (indeed an infinite series 
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of ambiguities) occurs. Note that, what is in fact gencrahy expected (and iUustrated in a 
few exphcit calculations in the (1+1) dimensional sigma model[8, 9]) is that the complete 
NP result is unambiguous, but the separation of its perturbative from its NP "operator 
product expansion" (OPE) contribution, is not. More precisely the two contributions are 
intrinsically related through the need of a definite prescription in choosing the integration 
path in the Borcl plane, if using e.g. Borel resummation methods, necessary to avoid the 
renormalons in both parts and to rcsum the otherwise ill-defined factorially divergent per- 
turbative series. Now in contrast, as we will show, a more interesting aspect of our results 
is that no such prescription is needed for the NP infrared renormalon disappearance in the 
GN mass gap. More precisely, the mechanism is such that the perturbative (re)expansion 
(to be specified) of the pole mass at order 1/A^ defines a directly Borcl summable new 
perturbative series, which can be thus unambiguously separated from the NP part. 

The paper is organized as follows: in section 2, briefly recalling some known basic 
results of the GN calculations at next-to-leading 1/N order, we shall define more precisely 
our "naive" perturbative expansion of the 1/N order GN pole mass and exhibit its infrared 
renormalon divergences. In order to examine those issues, we will consider first the GN 
model with a non-zero mass term, and single out a contribution which corresponds to a 
purely perturbative information. (The ordinary perturbative expansion of the true mass 
gap being of course zero to all orders, since the mass gap is non-perturbative and of order 
A ~ e"'^'^/^^^^ where g is the coupling). One can consistently recover the true mass gap in 
the massless Lagrangian limit. In section 3, we rederive the known [10, 4] exact calculation 
of the GN mass gap at 1/N order, but in the light of the results of section 2. Namely, we 
stress the difi'erences between this exact calculation, which automatically takes into 
account the full NP contributions, and the naive pole mass expansion of section 2 based 
on a purely perturbative information. We then perform in Section 4 a detailed calculation, 
using a convenient (Mellin transform) method, to exhibit the precise interplay mechanism 
between purely perturbative and NP contributions within the 1/N expansion, and leading 
to the expected absence of renormalons in the complete NP mass gap. We give in section 
5 conclusion and some more speculative comments on a possible generalization of those 
results to QCD. Finally, two short appendices are devoted to additional technical details. 

2. Perturbative renormalons in the 1/N naive pole mass 

We start by briefiy recalling the standard construction[10, 4] of the mass gap of the 0{2N) ^ 
GN model at order 1/A^, with a slightly different approach. Here we consider in fact the 
model with the addition of a Lagrangian mass term, in order to define the pole mass in a 
somewhat more universal manner, making in particular the link with analogous quantities 
in other models (QCD typically) more transparent. As usual the true mass gap is to be 
considered in the chiral, massless Lagrangian limit. The Lagrangian thus reads 

Cgn ='^i^^ -m^^ + (2.1) 

^From now on, all expressions correspond to the 0{2N) model, for easiest comparison with the exact[4] 
and 1/Af[6] mass gap results. 
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and the graphs contributing to the two-point function at next-to-leading order are 
shown in Fig. 1. At this 1/N order, the renormahzation procedure is relatively simple, since 
only the mass term is renormalized, which can be most simply performed by subtraction 
of the most divergent terms in a taylor expansion of the integrands [10]. The expression 



Figure 1: The GN mass graphs at order 1/A''. 
corresponding to Fig.l reads, after angular integration: 
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with M = A = lie ^""/(a^sCm)) i\^q mass gap at leading 1/N order (/x is the arbitrary 
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with A = {q^ + + M|,)2 - Ap^ M|, {q is Euchdean but p Minkowskian) , ( = {1 + 
A.Mp/q^Y^'^ , and the dressed propagator (wavy line) is 



G[q'] = 



gN.. X + i. 



+ 0(l/iV) . (2.4) 



To obtain the pole mass at the ncxt-to-lcading order, and in the masslcss limit, it is 
sufficient to replace ^ = Mp = M[l + 0{l/N)] in Eq.(2.3),(2.4) up to l/iV^ ^^^^^^ This 
reads 
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where a "factorization" scale n regularizing the integrals was introduced. The renormahza- 
tion by subtraction will remove the divergent part of integrands, also regularized in term 
of /X. This sharp momentum cut-off procedure is rather similar to the one in ref. [6] , giving 
results perfectly consistent with these authors, as we will see in next section, and is conve- 
nient for our purpose here, where we are primarily interested in the asymptotic behaviour 
of (2.5) (thus for /Lt 3> A, but kept finite). 
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For simplicity wc let aside in this section the tadpole second graph in Fig. 1, as its 
exact expression in (2.5) does not give any factorially divergent perturbative coefficients 
when expanded^. [It gives however a finite contribution to the mass gap, as we shall see in 
next section when a precise evaluation of the mass gap and its asymptotic behaviour will 
be considered.] We shall next examine in some detail how a naive perturbative expansion 
of (2.5) leads to infrared renormalons and related 0(A//Lt) ambiguities. Before to proceed 
to see that, it is perhaps important to better explain what we exactly mean by "naive" 
perturbative expansion of (2.5). Clearly, the expression (2.5) is exact and complete at 
1/7V order, which means implicitly that it contains a priori both perturbative and NP 
contributions. One may argue that, in fact, the ordinary perturbative expansion of the 
true mass gap is zero to all orders, the mass gap being a NP quantity of order A. However, 
this is true strictly speaking only when considering the massless GN model, i.e. for m ^ 
in Eq. (2.1). Alternatively, one may consider the calculation at this order of Eq. (2.5) 
as valid for the massive GN Lagrangian case where say, m ~ Mp = 0(A) (for completeness, 
the case where m S> A is also briefly discussed in Appendix A). Remark indeed that in 
(2.5), the 1 — C factor, which corresponds to the skeletal "one-loop" integrand of the 
first graph in Fig. 1, is essentially obtained from the pole mass condition p = Mp, 
irrespectively of whether Mp is 0{A) or not. The remaining (^-dependence is the dressed 
scalar propagator, as defined in Eq. (2.4), where the leading order mass gap definition 
is only used in the right-hand side of Eq. (2.4). Thus the perturbative form of (2.6) is 
quite generic, formally also applying to e.g. the QCD pole quark mass [but with obviously 
a different g^-dependence replacing the 1 — C term] . Now, since in the present case one can 
control everything at the order explicitly, our purpose is to consider the kind of results 
obtained when the NP knowledge is partially omitted. More precisely, let us assume that 
within the next-to- leading 1/N graph structure in Fig. 1, only the purely perturbative 
part of the dressed scalar propagator is known. Our motivation is that this is analogous 
to the situation in a more complicated theory like QCD typically. Starting thus from the 
"purely perturbative" information means replacing in (2.5) the dressed scalar propagator 
by the effective coupling^: 



where in the bracket of Eq. (2.6) the 0{M'^ /q^) "NP" corrections to the effective coupling 
are of course exactly known in the GN case and can be obtained explicitly here from a 
systematic expansion of the right-hand side of Eq. (2.4) for large q^. Again, we stress that 
omitting this NP part is evidently not fully consistent, but the purpose of this (somewhat 
artificial) separation of the different contributions is to examine precisely how the renor- 
malons (would) appear, in a way very similar to the QCD pole mass, and more importantly 

^Due to the second equality in (2.4), valid at this order, which makes the tadpole graph simplifying 
to purely polynomial integrals. Beyond the 1/N order, one could still in principle choose the arbitrary scale 
/Lt such that (2.4) holds. 

^In the sequel we rescale hog = N/{2n) g g, to define the 1/N expansion properly, and absorbing as 
well the 2n factor just for convenience. 




(2.6) 
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from which detailed mechanism do they finally disappear in the GN case, once the complete 
NP information is taken into account. 

The simplest standard procedure to exhibit the IR renormalons is by expanding the 1 — C 
term for small 

and expanding the effective coupling ^(g^) at one loop order of renormalization group (RG) 
in powers of bog{iJ,) ln[(/^//x^]. From (2.7) it is seen that the leading singularity comes from 
the term which, combined with the (^(/x) ln[q'^//x^])" terms of the effective coupling 
expansion, produces factorially divergent perturbative coefficients: 
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The non sign- alternation of those factorial coefficients implies [7] that the corresponding 
series is not Borel summable: the Borel integral corresponding to (2.8) reads 



(2.9) 



so that the pole at to = on the integration path gives the leading ambiguity for the 
pole mass 



/oo-|-'ie i^oo—te 
- ]{dt e-'/\l - 2t)-^) = ±i7rMe-*°/5 (x A . (2.10) 

Note that this leading renormalon ambiguity of C(A) is completely similar to the one de- 
rived for the quark pole mass in QCD[5]. The fact that the perturbative contribution Eq. 
(2.9) is proportional to the factorization scale fi may appear at first puzzling (since it is 
supposed to be a contribution to the mass gap ~ A ^ /i) , it is in fact also a known feature 
of the non Borel summability of the corresponding series[7] with this explicit cutoff regu- 
larization: namely, it indicates us that an additional NP contribution, also proportionnal 
to the factorization scale, is necessarily needed (at least in this regular ization scheme) in 
order to recover the right mass gap dependence ~ A. 

Although our above derivation (starting from purely perturbative information) is very stan- 
dard, the fact that the GN model perturbative pole mass at order also has the specific 
structure (2.8)-(2.10) of infrared renormalons, was perhaps not clearly appreciated before, 
to our knowledge^. Similarly, we can easily check from (2.3) that the leading asymptotic 
behaviour of perturbative coefficients of the two-point function for S> is ~ X^pp! g^-, 
which accordingly gives a less severe ambiguity of 0{K^ /p^), where again one can note the 
similarity with the QCD off-shell p^ ^ quark correlation function case [5, 7]. 
It is in fact possible to go a step further and to characterize at arbitrary next orders 



*In ref. [1] appeared already the (earliest) discussion on ultraviolet renormalons in a field theory frame- 
work, which were shown to be harmless (Borel summable) for the asymptotically free GN model. 
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the renormalon properties of the naive perturbative expansion of the GN mass gap. 
Consider the second order (exact) RG dependence of the effective couphng[ll]: 
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(2.11) 



defining / recursively, with g = g{n) and the beta function f3{g) = — 2605^ — 26i_g'' — • • • 
(where for clarity we reintroduce the original coupling and RG coefficients, i.e. before 
rescaling of g). We can put Eq. (2.5), with (2.11), directly into the form of a Borel integral, 
after a convenient change of variable[12], defining the (Borel) variable t as 
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Taking expression (2.6), but using now Eq.(2.11), we find after some algebra 
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+ 2 2^p!r[3/2-p] 



where C = 6i/(26o) and a = 1/5 + 61/60. To obtain (2.13) we expanded Eq.(2.7) in powers 
of /M'^, and used Eqs. (2.11,2.12). This gives the complete (leading and all subleading 
orders) series of infrared renormalons (initially corresponding to cuts at botp = 1/2, 1, ..{2p+ 
l)/2, p G N*). The change of variable (2.12) makes calculations more convenient since 
expression (2.13) has only a cut at t > I/60. Now we can calculate the ambiguity to 
all (perturbative) orders (of course still limited to 1/N order), which we define by the 
difference of contour above and below the cut. We find, again after some algebra: 
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and identified the MS scale A = /le 



-l/(2M)f 



[bog] "^[1 + (61/60)5]*^ consistently at second 



RG order. Eq.(2.14) thus gives the full series of ambiguities due to infrared renormalons 
(for the first graph in Fig. 1), in the form of power corrections of order A^, with the first 
term in the bracket of (2.14) the leading order ambiguity of 0{A). If 61 = (e.g. at first 
RG order) expression inside the bracket of Eq. (2.10) simplifies to (1 + r^/4)^/^ — r/2, 
r = A/M (indeed sufficient at this 1/N order, since from Eq. (2.2) ri is already the 1/N 
term.) 
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As already emphasized, the above derivation of (2.14) only uses information that is in fact 
purely perturbative: the effective coupling at second RG order, Eq. (2.11), and the specific 
GN mass kinematic dependence, Eq. (2.7) of the (one-loop) skeletal first graph in Fig. 1. 
Accordingly, a similar derivation is possible for the QCD quark pole mass renormalon prop- 
erties (indeed the equivalent of the information in (2.14) in the QCD case is also known, 
though perhaps expressed in a slightly different form[7]). 

Actually, by considering only the first graph of Fig. 1 we slightly oversimplified the com- 
plete renormalon picture for the GN model: clearly, the second tadpole graph also gives 
renormalons, if considered purely pcrturbatively. These are easily analyzed similarly to 
Eqs. (2.6)-(2.7), and the expanded integrand for q^ ^ leads to renormalon poles at 
^o^p = p G N*. Incidentally, the pole at hot = 1 exactly cancels the one in (2.13), but the 
leading as well as all subleading poles for hotp > 3/2 in (2.13) remain uncancelled. 
In summary, we thus observe from the structure of (2.14) that even resumming the full 
integrand 1 — C (taking the full series of sub-leading renormalons) does not remove in 
any way the leading ambiguities (even if there are some "accidental" cancellations among 
subleading poles at this purely perturbative level, between the two graphs of Fig. 1 as 
above discussed). This is not at all surprising, since as emphasized (2. 13), (2. 14) are still 
perturbative calculations. This ambiguity in the purely perturbative piece of the pole mass 
is severe, being of the same 0{A) order than the mass itself, similarly to the QCD pole 
mass[5] ambiguity. As explained above in detail, this originates from the 1 — C factor, spe- 
cific to the pole mass, which has an expansion for — whose leading term is 0(1^1"^), 
cf. Eq. (2.7). It is interesting to compare those pole mass results with what happens for 
other related quantities of interest. For instance, the GN fermion condensate {^^) (which 
is related to the mass gap, since e.g. {^^) = NMp in the N ^ oo limit), has a leading 
ambiguity in its perturbative expansion which is much power suppressed. More precisely, 
the relevant part of the condensate can be obtained, roughly speaking, by closing the ex- 
ternal leg of the first graph in Fig. 1. A straightforward calculation along the lines of Eqs. 
(2.6)-(2.10) gives for the relevant integral: 



where cq and ci are unessential constants not specified here. Accordingly, when takin 
in (2.16) only the perturbative part g{q^) of the propagator G[q^], cf. Eq. (2.6), one 
finds for the Borel transform an IR renormalon pole at hot = 1, thus associated with an 
ambiguity of 0(A//x)^. This result for the IR renormalon in the perturbative tail of the 
GN condensate was in fact obtained long ago[13], both in this heuristic form and rigorously 
by using a lattice regularization with Wilson fermions, where it was precisely shown that 
this absence of a leading order renormalon is due to the chiral symmetry. As concluded by 
these authors, even if the perturbative expansion of the condensate is not, strictly speaking, 
Borel summable, the ambiguity may be neglected since it is power suppressed with respect 
to the NP condensate itself, the latter being of 0{A). A similar suppression due to chiral 
symmetry of the leading IR renormalon also occurs for the QCD quark condensate in four 
dimensions. 




(2.16) 
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Coming back to the GN pole mass case, we have however to stress once more at this point 
that all the perturbative renormalon results derived in this section are in a sense artifacts, 
since the perturbative expansion that we considered was on purpose artificially generated 
as an incomplete part of the full l/N expression of the mass gap. 

3. Borel summability of the exact 0{1/N) mass 

Alternatively, since expression (2.5) is exact at 1/N order, we can calculate exactly the 
expression of the mass gap, i.e. without truncating (2.5) to its perturbative expansion. 
Indeed integral (2.5) can be evaluated analytically exactly: after a convenient change of 
variable = tanh(|) we obtain 



withx = (1 + 4MV/x2)V2 = i/tanh(0/2) (i.e. = ln[(x + l)/(x - 1)] > 0), and Ei(-x) = 

— dte"^ /t {x > 0) the Exponential Integral function. The term — 21n(cosh[^/2]) in (3.1) 
corresponds to the (unsubtracted) tadpole graph of Fig.l, and the terms lnln(/x^/M^) = 

— Ing and ln//^/M^ = 1/5 are the subtraction terms for the first and tadpole graphs, re- 
spectively. One can easily check the finiteness of (3.1), if letting the "cutoff' /i oo. Now, 
as already stressed, here we are interested in the complete asymptotic behaviour, thus let- 
ting 3> A but kept finite, retaining eventually all the power correction terms in A///. The 
function Ei{—x) for a; > has an asymptotic expansion with factorial but sign alternating 
coefficients, therefore explicitly Borel summable and perturbatively unambiguous. More 
precisely when re-expanding (3.1) in perturbation, we obtain 



where the higher order power correction terms, that we do not give explicitly here, can be 
obtained by a systematic expansion in M^///^ of 9. The tadpole graph in Fig.l contributes 
a constant term 2 In 2, relevant to the mass gap determination, but does not contribute the 
factorial asymptotic behaviour of the perturbative series, as already mentioned. 
We emphasize that the specific 1 — C form of the integrand in (2.5) plays an essential role 
for the Borel summability of expression (3.2), which is accordingly peculiar to the pole 
mass. In contrast, the off-shell two-point function expression (2.3) for arbitrary may 
be evaluated similarly non-perturbatively (at 1/A^ order), and does not lead to a directly 
Borel summable perturbative series, in consistency with the results obtained and discussed 
previously in ref. [10]^. We shall come back to this specificity of the pole mass in more 
details in next section 4. 

We also stress that the use of the 0{1/N) mass gap relation Eq.(2.2) in the Lagrangian 
mass m — > limit, though it greatly simplifies the analytic evaluation of Eq.(2.5), plays no 

^ Technically, the simplest such terms can be expressed in terms of Ei{d), which accordingly (since 
9 > 0) has asymptotic expansion with same sign factorial coefficients and imaginary part 
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particular role in the good asymptotic properties of Eq.(3.2): more precisely, starting from 
the exact 1/N mass expression (2.5), and introducing the arbitrary mass dependence m, 
Borel summability is maintained with an asymptotic expansion similar to Eq.(3.2) for any 
value of the pole mass Mp S> A (see Appendix A). 

Finally to obtain the correct mass gap M^/A at next-to-leading 1/N order from expression 

(3.1) , one should yet introduce the MS basic scale above defined after Eq. (2.15). Dropping 
terms of higher 0{1/N'^) order, we obtain 

, ,p T^ 1 N , l-|-21n2 — 7e; , , 

M^/A = 1 + ^ + O(^) = 1 + (3.3) 

in agreement with ref.[6]. Note however that our expression (3.1) differs in fact from rcf.[6], 
more precisely by the Ei[—6] term. This is simply because in ref.[6] all terms vanishing as 
inverse powers of ^ were dropped, which is sufficient to identify the mass gap Eq. (3.3). 
From our result, this is consistent because at the (non-perturbative) level of Eq.(3.1), all 
those "power correction" contributions from Ei{—9) can be unambiguously separated, thus 
dropped from Eq.(3.2), to let only the part relevant to determine the mass gap. The explicit 
Borel summability of the "educated" perturbative expansion of the pole mass at 1 /N order, 
Eq. (3.2), confirms the consistency of the whole procedure. An equivalent signature of Borel 
summability here is the absence of positive powers of in (3.2), in contrast with the purely 
perturbative contribution Eq. (2.8). Applying the same procedure to other quantities than 
the pole mass, one eventually ends up with asymptotic expansions with non sign-alternated 
factorial coefHcients and /x^ power terms with P > (as is clear from the above expression 
of 9, cf. footnote 5). Within a complete NP calculation, the latter fi^ terms in fact cancel 
when combining the OPE and perturbative parts, as illustrated in explicit calculations for 
some vacuum expectation value[8, 14] and off-shell correlation functions[9] of the 0{N) 
sigma model at 1/N order. As mentioned in introduction, this also illustrates that, though 
any complete NP calculation is expected to be renormalon ambiguity free, one cannot 
expect in general (for an arbitrary Green function) to be able to single out unambiguously 
the perturbative contributions. 

Coming back to our result Eq. (3.2), it appears thus immediately in apparent contradiction 
with those obtained starting from purely pertubative expansions, (2.8-2.14) above. We 
shall examine in next section how to reconcile these two different pictures. 

4. Explicit disappeairance of IR renormalons 

How exactly the "bad" factorial coefficients with no sign alternation in Eq. (2.8) disappear, 
or more precisely transmute into "good" sign-alternated factorials, in the exact expression 

(3.2) ? Clearly, it can only be that this necessary change of the wrong sign perturbative 
expansion factorial coefficients should occur with the NP power expansion contributions: 
the weak point of the standard perturbative renormalon picture is that we have expanded 
the 1 — C term in the infrared low q'^ regime, while keeping the short distance, perturbative 
effective coupling form of the propagator. This is motivated, as already mentioned, from 
the fact that the latter information is a priori the only accessible one in more involved 



-9- 



theories such as QCD. In the present case, as we know exactly the mass at order one 
may at first hope to infer such transformations by examining e.g. the systematic short 
distance — oo and long distance q"^ —>■ power expansion of the integrand in (2.6), 
which are perfectly well-defined^. 

In fact, to see the vanishing of renormalons operating needs a little more sophisticated 
analysis. Following e.g. refs.[8, 9], we introduce the Mellin-Barnes (MB) transform for a 
part of the integral Eq. (2.5), which then takes the form (again omitting here the tadpole 
graph for simpler illustration): 



where the Kernel (inverse MB transform) is defined in our case by 



(4.1) 



K(s,t)= r 



C-1 

c+i 



(4.2) 



The MB transform method main purpose is that it will exhibit precisely the singularities 
of the integrand, in the Borel plane t of interest. The sequel is just algebraic manipulation. 
Changing again variable C,~^ = tanh(^), except for the 1 — C term kept on purpose, using 
Eq. (2.7), as an expansion in g/M, (4.2) can be evaluated exactly, and 
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-a r[l + a-2g]r[-a/2 + g + t] , 

" T[l + a/2-s + t]{l + s) ^ ' 



defined for i?e[s + t] > a/2 and i?e[2s] < 1 + a. The latter conditions are such that integral 

(4.2) converges, and play essential role in determining the singularities. The variable a in 

(4.3) is simply the power oiq/M ~ sinh0 in expansion (2.7), with coefficient Ca respectively. 
Thus a = — 1 with c_i = —2 corresponds to the leading renormalon, and a = 0, ..2p + 1, 
p G N to subleading ones. To evaluate the s integral one can simply close the contour 
on the left, and sum over residues of the poles included in this domain (since = 
(M^/f/^)^* decays exponentially fast for the asymptotic regime q^ ^> we arc interested 
in). Expression (4.3) has (simple) poles at s = —1, s = a/2 — t — k, and 2s = a + 1 + k, 
A; G N, where the latter poles do not contribute for the relevant contour. The final result 
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where the first term in bracket corresponds to the residue of the pole at s = —1, while 
other terms correspond to the sum over residues of the poles at s = a/2 — t — k, and 



®E.g. the infrared 
dependence: [Clii[|^]]~^ 



expansion of the dressed propagator is purely polynomial, with no log. 
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we used (M^//i^)* = e~*/^. In (4.4) one sees that the first term in the bracket and the 
summed terms both have poles at t{a, fc) = 1 + a/2 — k, which can occur at i > depend- 
ing on a, k values. On more physical grounds, the contributions from the first term, the 
initially s = —1 pole, originate from power terms (M^/g^)~* in Eq.(4.1), and correspond 
intuitively to non-perturbative "OPE" contributions, while the perturbative contributions 
are those multiplied by e~*/^ in (4.4). Even though in fact both contributions are in the 
end of the same (non-perturbative) order, such as to compensate the bad poles in t > 0, 
it is natural within the Borel transform formalism to distinguish those two contributions 
as "non-perturbative" and "perturbative" respectively, since typically the second contri- 
bution has the characteristic t-dependent overall factor e~*/^ = (M^//x^)* multiplying an 
expression with poles at t > 0, so that when this latter contribution is considered alone it 
organizes as an ordinary perturbative expansion when integrated over t. More precisely, 
for example the poles at t = 1/2 for a = —1 correspond to the leading order renormalon, 
with k = 0. Indeed, keeping only the leading renormalon perturbative terms, oc e~*/^, for 
k = 0, one recovers exactly Eq. (2.9). 

Now it can be easily checked that this t = 1/2 pole in fact cancels exactly against the first 
NP term t = 1/2 pole, and similarly for all sublcading poles at t = 1, 3/2, ..(2p+ l)/2. This 
is the announced disappearance of renormalons. Moreover all cancellations happen, for a 
given pole at t = 1 + a/2 — k, between NP and perturbative terms of the same k values. 
Since (4.4) is in the form of a Borel integral, and after cancellations all remaining poles 
occur at t <0, it defines a Borel summable series, whose leading terms just correspond to 
the asymptotic series defined in (3.2). To see that, it is simpler to alternatively proceed 
directly with (4.2) using the change of variable C, ^ (f) for any terms in the integrand, i.e. 
without going through power expansion (2.7): after MB transform one ends up directly 
with a Borel integral where no poles at all appear at i > 0, and which exactly gives the 
asymptotic series in Eq.(3.2), including correct finite terms —^e + ••• [NB there are also 
poles at t = in (4.4), which as usual[7] simply corresponds to the UV divergences, and 
are removed consistently by the appropriate subtraction terms, that we do not display 
explicitly.] For instance, illustrating only the terms from the first graph in Fig.l, after 
cancellations of the t > poles the MB transformation gives 



^2 f'^ 



M2 1 ^ 1M2 3 + 2t m2 ,^,M^o. 



1 /tf2 °° 

~-^[7^ + -^5E(-l)W] (4.5) 

n=0 

to compare with Eq. (3.2) (where we used e.g. dtt^^[l/{t + l) — e^^] =^e)- Note that in 
Eq. (4.5) all the "bad" positive powers of jJi/M, that signalled the non-Borel summability in 
the corresponding incomplete purely perturbative expressions Eqs. (2.9), (2.13), have now 
disappear. Of course, proceeding in this "direct" way is nothing but a consistency check 
that the MB transform gives a correct alternative calculation of the asymptotic expansion 
of the exact integral Eq.(2.5), which wc started from anyway. But our explicit separation 
of the expanded perturbative renormalon part (2.7) in connection with the MB method 
allows to visualize explicitly the renormalon disappearance order by order in (4.4). For 
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completeness note that a very similar MB transform analysis can be performed for the 
tadpole graph renormalons, which disappear similarly from a detailed analysis that we do 
not display here. 

It may be of interest to examine the influence of the regularization scheme on the 
results derived above. As expected, the renormalon properties are independent of the reg- 
ularization used, but they manifest themselves in a quite different form e.g. in dimensional 
regularization [8, 7]. The momentum cutoff regularization prescription considered here, in- 
troducing an explicit factorization scale /i, cf. Eq. (2.5), is clearly very convenient for our 
purpose since it directly displays the NP power correction terms, while the interpretation 
of the latter for instance in dimensional regularization (DR) is not completelly obvious. 
Nevertheless, the appearance (and disappearance) of IR renormalons as discussed here can 
also be seen equivalently in DR. As is known from general analysis[8, 7], the IR renor- 
malon ambiguities in DR manisfest themselves in the fact that the Borel integral over the 
perturbative contributions has to be evaluated for complex e, and gives a different result 
depending on whether /m[e] > or /m[e] < 0. The difference between the two results 
is proportional to the NP contribution. In a rather simplified analysis of the GN pole 
mass (see Appendix B for more detail), considering the equivalent of the MB plus Borel 
transform Eqs. (4.1)-(4.4) in DR, one arrives at an expression quite similar to Eq. (4.4), 
but in which now both terms have poles at t = a/2 -I- 1 — — e, and their residue cancel 
each others order by order in a way completely equivalent to what was discussed above for 

To conclude this section, going back to the more convenient cutoff regularization, one 
may alternatively understand perhaps more qualitatively the specific pole mass renormalon 
disappearance by examining the asymptotic behaviour of the off-shell mass expression (2.3), 
first expanded in powers of M^/p^, and proceeding with the MB transform similarly to 
the above described procedure. [This is then a completelly similar calculation than the 
one performed in details in ref.[9] for the 0{N) sigma model.] Skipping many details, we 
simply indicate sketchily that for the GN model, it would give coefficients of the (M^/p^)" 
terms with a form similar to Eq. (4.3), but with essentially the replacements P^', 
(1 -|- s)~^ Sj(l + s -|- < i < n — 1. Those poles at s = — 1 — i, of residues 
~ r[t — 1 — i], produce in turn poles at (integer) t = 1,2, ..n. Again, the integral over t 
may be cast into a Borel transform formally similar to (4.4), identifying a NP part and a 
perturbative part. However, in this off-shell case, the specific poles at to = 1, 2, .. within 
perturbative terms of order (M^/p^)", are cancelled by NP term of different order n + to. 
Since individual terms are singular, one cannot truncate the power expansion unless a 
definite integration contour prescription to avoid the poles is defined [9]. In other words, 
the separation between the NP and perturbative part in this / case is ambiguous, 
which means the non Borel summability of the perturbative series of the general two-point 
function (2.3) for arbitrary]?^, as obtained from a direct calculation, cf. remarks in section 
3 (footnote 5). In contrast, for the pole mass any (M^/p^)" terms are replaced by 1, which 
"fiatten" all orders of the M'^/p^ expansion, so that the different cancellations of t > 
poles now occur all at once. 
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5. Conclusion and perspectives for QCD? 

In this paper we have exhibited in details the non-trivial disappearance of the "naive" 
perturbative IR renormalons of the GN mass gap at order implying the direct Borel 
summability of the "educated" perturbative expansion, as defined by Eq.(3.2). Given 
the detailed rcnormalon vanishing mechanism, we are also confident that it should work 
similarly beyond order, though an explicit check has not been attempted. Concerning 
strictly our 1/N results for the GN pole mass, a number of additional remarks may be 
useful. 

-First, we should perhaps again emphasize that the result Eq.(3.2) may be not a surprising 
one, as it could have been easily extracted from previous analysis of e.g. refs. [10, 6], if 
not explicitly displayed there. As motivated in introduction, our main purpose here was 
to illustrate in a calculable model how the NP contributions to the pole mass organize 
to eliminate completely the renormalon artifacts, even though the latter are somewhat 
artificially introduced by defining the naive perturbation theory, similar to the perturbative 
expansion of QCD quark pole masses. 

-Second, the (would-be) perturbative ambiguities due to IR renormalons for the pole mass 
here discussed are severe, being of the same {0{K) order than the pole mass itself. In 
contrast, there are other related quantities of interest, like most typically the GN model 
fermion condensate, where the leading order potentially dangerous renormalon singularities 
in the perturbative part disappear[13] by explicit cancellations due to the chiral symmetry, 
as reexamined briefly in the end of section 3. There, the cancellation of renormalons is such 
that the actual perturbative ambiguities are of higher order than the condensate dimension, 
thus power suppressed and to be considered harmless. Though the latter results are not 
unrelated with the present study, note that the IR renormalon cancellation mechanism as 
we exhibited here is of a rather different nature, since all (but not only the leading) IR 
renormalon singularities are annihilated each other between the purely perturbative and 
NP contributions, see e.g. Eq. (4.4). But we stress again that these results are very specific 
to the pole mass. 

It is thus tempting to speculate briefly on the possibility of ultimate (NP) renormalon 
disappearance in QCD similarly (to some extent) to the mechanism discussed here. QCD 
in the massless quark limit also has a mass gap, since the approximate chiral symmetry 
of the light quark sector is dynamically broken. Now if we assume as usual[7] that the 
dominant IR renormalon contributions to the pole mass is given by graphs of a topology 
analogous to the first one in Fig 1, but with a gluon propagator (wavy line) dressed now 
with massive, constituent quarks of mass Mq ^ Aqcd = A, one may expect this propaga- 
tor to behave in the infrared in a way similar to Eq.(2.4). (Assuming also that such quark 
loops are complemented with appropriate QCD gauge sector contributions, so to match 
the correct beta function in a gauge invariant way, as indeed usually assumed in most 
QCD infrared renormalon issues[7]). We see no reason why such assumption on the NP 
behaviour would not be consistent with the usual perturbative behaviour, and in particular 
with the standard heavy quark pole mass renormalon picture[5], for which Mq ^ A. This 
is also irrespective of the fact that confinement in QCD ultimately makes the pole quark 
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mass relevance somewhat elusive. What can be still theoretically relevant would be to have 
in this way a procedure to evaluate the NP contributions to the light constituent quark 
masses Mq ~ A from first principles (and perhaps more interestingly for the NP order 
parameters related to chiral symmetry breaking). A detailed QCD analysis is however 
beyond the present scope and left for future work. 

Let us finally mention that, indeed, there exist more "direct" ways of modifying the asymp- 
totic properties of the perturbative expansion of e.g. the mass gap, generically in asymp- 
totically free models. It relies only on the purely perturbative information, but is based on 
a modification of the usual perturbative series. Such a method[15] appears to bypass the 
explicit disappearance of renormalons here exhibited, by directly removing the perturba- 
tive renormalon divergences, at least for adequate ranges of an (alternative) perturbative 
expansion parameter values compatible with the chiral limit. 



A. Asymptotic behaviour for Mp ^ A 

As mentioned in section 3 the use of the mass gap relation Eq.(2.2), strictly valid only at 
0{1/N) and in the chiral limit (Lagrangian mass m ^ 0), plays in fact no particular role 
in the asymptotic properties Eq.(3.2), of the exact 1/N pole mass. Consider Eq.(2.5), but 
now for Mp 3> A (which illustrates e.g. the massive theory case with Mp ~ m 3> A): it 
takes the form 

n = / #[e-« cosh I - |(1 + ,n(^))-'l(ln(^) + -A^)"' (A.l) 

to be expanded in powers of ln~'^(Mp/A^), where we used again (j) = ln[((^ -|- 1)/(C — !)]• 
Each coefficient of such an expansion contains a logarithmic divergence, renormalizable by 
subtraction (removing essentially the divergent (l-|-^)/2 piece in the bracket of Eq. (A.l)). 
After some algebra we find finally that the renormalized series has the leading asymptotic 
behaviour 

oo I 

rf^{Mp » A) ~ ^(-l)f+i^_^ (A.2) 

which qualitatively agrees asymptotically with Eq. (3.2), provided one identifies Mp ~ 
/X S> A, i.e. ln(Mp/A^) ~ 1/g. The sign alternation in (A.2), leading to Borel summabil- 
ity, again makes the main difference with the behaviour obtained starting from the naive 
perturbative analysis of section 3. 



B. IR renormalons (dis) appearance in dimensional regular izat ion 

In dimensional regularization (DR) with D = 2 — e (so that the GN model is renormalizable 
for Re[e] > 0), Eq. (2.5) may be cast in the form (omitting once more the tadpole graph 
to simplify): 
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where ^ = (1 + 4M^/g^)(^/^\ and the DR regularized auxihary field propagator 
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(l-e)r[e/2] / dx{M^ + q^x{l 
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+ 0{1/N) (B.2) 



has the e — limit Eq. (2.4) (of course after coupling renormalization go = Zgg, which 
details are not explicitly displayed, absorbing the divergence in r[e/2]). Skipping the 
detailed calculation and unessential overall factors, the important point is that the large 
q'^ behaviour of Eq. (B.2) in DR, equivalent of Eq. (2.6), replaces integrals over Inq^ with 
integral of the form^ 



a>-l I 



/ (F' xa/2-c/2-c//2 



(B.3) 



(where the sum over a comes from the expansion of 1 — C, similarly to Eq. (2.7), while the 
sum over I simply comes from the geometic expansion of Eq. (B.2)). Though integration 
of (B.3) for finite e does not lead to factorial divergences unlike its e — >^ equivalent Eq. 
(2.6), it gives poles at 1 + a/2 — e(l + l)/2 to be evaluated by taking e complex, and the 
e ^ result diff'ers for /m[e] > or /m[e] < 0. This is similar to general results[8, 7] on the 
way in which renormalon ambiguities manifest themselves in DR. Going now to Mellin and 
Borel transform as examined for the cutoff regularization in section 4, the DR equivalent 
of Eq. (4.4) takes the form ^ 
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(B.4) 



where the first and second term in bracket occur as residues of poles at s = — 1 + e/2 and 
s = a/2 — t — k, respectively, where s is the Mellin transform variable, cf. Eq. (4.1). Both 
terms have poles at i = a/2 + l — fc — e/2, and their residue cancel each others order by 
order in a way completely similar to Eq. (4.4) for e ^ 0. 
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